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C^ , Abstract. The quasiparticle excitation spectra of a Bose gas trapped in a highly 

anisotropic trap are studied with respect to varying total number of particles by 
numerically solving an effective one-dimensional (ID) Gross-Pitaevskii (GP) equation 
C^ I proposed by Meteo et al. [1] recently. We get static properties and Bogoliubov spectra 

of the system in a large energy domain. We show that this method is computational 

' , . efheient and very acute for a condensate system undergoing ID to three-dimensional 

^ I (3D) cigar-shaped phase transition by comparing our results with those calculated by 

O ' the 3D-GP equation and analytical results obtained under some limiting cases. We 

I ^1 , identify the applicable parameter space for the effective ID-GP equation and find that 

this equation fails to describe a system with large number of atoms. We also identify 
that the description of the transition from ID BEG to 3D cigar-shaped BEG using 
this equation is not smooth as mentioned in Ref. [1] , which highlight the fact that for a 
C^ ' finite value of a±/as the junction between the ID and 3D limiting cases is not perfect 
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1. Introduction 

The character of elementary excitations is very important for understanding the 
macroscopic quantum behavior of a trapped Bose-condensed gas. Measurement of 
the collective modes in trapped gases of alkali- metal atoms [3] were carried out soon 
after the discovery of Bose- Einstein condensates (BECs), followed by various theoretical 
investigations and numerical analysis[4, 5, 6, 7, 8, 9, 10, 11, 12], especially about the 
low-energy excitations which are essentially of collective character, with respect to 
different traps. The full tunability of magnetic and optical traps opens an extraordinary 
opportunity to study in practice not only one- and two-dimensional (ID, 2D) Bose 
systems, but also dimensional crossovers influenced by the number of particles, size and 
shape of the system, interaction strength, and temperature [13]. Particularly, in a very 
elongated harmonic trapping potential, the condensate undergos a transition from quasi- 
ID BEC to three-dimensional (3D) cigar-shaped BEC with an increase of the number 
of atoms, or the aspect ratio (anisotropy) of the trap, i.e. A = uj-^juj^, where uj^ and 
bj^ are the axial and the radial frequencies of the trap potential respectively. Rapid 
progress in experimental techniques have made it possible to increase the aspect ratio 
of the trap, up to 1/700 [14, 15], which makes many conflgurations possible. For many 
theoretical works and numerical simulations a general and accessible methodology for 
determining the static solutions and excitation frequencies of trapped BECs in deferent 
regimes is required. 

The Gross-Pitaevskii (GP) mean-fleld theory[16, 17] has proven to be an 
indispensable tool in both analyzing and predicting the outcome of experiments of dilute 
condensates in zero temperature limit, the 3D form of which is 



d^ij.t) r 927/;(r,t) 



VtraA^) \ gN\i,\^ V^(r,t), (1) 



dt 2m dx^ 

where g = A-nfi^as/m is the 3D coupling constant with a^ being the bulk s-wave 
scattering length and m being the mass of atoms. The 3D order parameter ijj{r,t) 
for a gas of N atoms is normalized to unity. Providing a highly anisotropic harmonic 
trap potential, which is separable with respect to the x, y and z axes 

Vtrap{r) = V±{x, y) + V{z) = -mu\x^ + -mw^y^ + -mulz^, (2) 

is applied, the motion of atoms along the strongly conflned directions is frozen and thus 
the degrees of freedom of the condensate in these dimensions is reduced. The condensate 
is in the ground state corresponding to the motion in these directions, since they do not 
have enough energy to reach the related excited states. If the axial frequency of the trap 
is much smaller than the radial frequency, i.e. Uz -C cuj^, the corresponding dynamics 
of the systems becomes effectively ID. The local density of the axial condensate, piD-, 
takes a one dimensional form by integrating over the transverse coordinates, 

p^D{z,t) = N j dh^\,p{v,t)\\ (3) 

and the condensate wave function can be factorized in the form 

V'(r, t) = 0x(rx, Pid{z, t))^{z, t), (4) 
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and pir, = 7V|0(2;,t)p is obtained providing the radial wave function 0_|_ is normalized 
to unity. 

In the ID limiting case the condensate dynamics is governed by the low dimensional 
GP equation where the coupling constant is gio = g/2iia\ = 2hu±as with aj_ = 
Jh/mu_L. This ID form GP equation is only applicable when the condition a^piD <^ 1 is 
satisfied, where the mean-field interaction energy can be treated as a weak perturbation 
and then the condensate wave function that minimizes the energy functional is given, 
to the lowest order, by the Gaussian ground state of the harmonic oscillator, and the 
condensate is tightly confined in the radial direction. This requirement is equivalent to 
a limitation on the condensate occupation number to A^ < aj_/as when huj± exceeds the 
magnitude of the mean-field interaction energy. 

The Gaussian approximation of the radial density profile becomes, however, 
poor when the number of particles or the aspect ratio of the trap increases, where 
the condensate will be quasi-lD or 3D cigar-shaped. To find an efficient way to 
describe the axial dynamics of this kind of condensate with reduced dimensionality 
several authors have followed different routes to introduce more accurate analytical 
approximations which relate the radial and axial wave function with the ID density 
profile[18, 19, 20, 21, 22, 2, 1]. The effective-GP equation has been proposed recently 
by Mateo et al. [1] is the most accurate and simple. By deforming the interatomic 
interaction one can get an effective ID-GP equation [2, 1] 

indt(t>{z, t) = -— V20(2, t) + [V{z) + hujJl + AasN\(t)\'' ]0(z, t) (5) 
Im ^ 

which describe the two limiting cases and the crossover between them very well. The 

basic idea in this approach is using a local chemical potential /i^ = Tiu ^^/T^^~cQ)l£l to 

include the contribution of the radial degrees of freedom to the axial dynamics. This 

provides us a very convenient tool to investigate the excitations of elongated BECs in 

different regimes. 

In this work we give a computational efficient way to calculate Bogoliubov 

excitation spectra of trapped BECs from ID to 3D cigar-shaped, including the crossover, 

configurations. By the Bogoliubov equations numerically based on the effective-GP 

equation 5 we get static properties and quasiparticle excitation spectra with a range 

of increasing number of atoms in a large energy domain. We identify the applicable 

parameter space for the effective ID-GP equation and find that the transition from the 

ID BEG to the 3D cigar-shaped BEG is not as smooth as mentioned in Ref. [2] which 

highlight the fact that for a finite value of a±/as the junction between the two domains 

is not perfect [11, 12]. In the low energy excitation domain the simulation results are 

compared with the analytical results. In Sec. 2 we investigate the ground state properties 

of the axial condensate obtained by numerically solving the effective ID-GP equation 5. 

We found the results are very accurate comparing with those calculated by the 3D-GP 

equation. In Sec. 3 we give the Bogoliubov equations derived from the effective ID-GP 

equation and the numerical methods we used. Analytical results for the very low lying 

excitations based on the variational method, hydrodynamical methods and the sum rule 
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approach are discussed in Sec. 4. In Sec. 5 we will analyze the numerical results of the 
Bogoliubov spectra comparing with the analytical results given in Sec. 4. 

In the numerical investigations we use az = Jh/mujz as the unit of length and l/w^ 
as the unit of time. Hence we define dimensionless wavevectors k = kaz, times t = uj±t 
to minimize numerical instabilities and simplify our expressions. 



2. Ground state properties 

Under Thomas-Fermi (TF) approximation the chemical potential obtained from the 
time-independent effective ID-GP equation is 

A4(l + iAZ=). (6) 

where Z is the dimensionless TF radius of the axial condensate which is defined by the 
value of z at which the equilibrium density Pid{z) vanishes. The density profile is then 
found to be 
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with pio = for \z\ > Z. Then from the normalization condition 
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_ piD{z)dz, 



(7) 
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we obtain the relation which determines the axial half length 



L^^z)' + liVxz)^ = nVx- 
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For a giving A the only relevant parameter is x- An approximate solution of this equation 
is given by [2] 
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_ (15x)^/^ + 1/3 57x + 345 (3x)^/^ 
In the limiting case % ^ 1 we can get 



(10) 



(11) 



This is the same result as we use 3D-TF approximation for a stationary 3D-GP equation, 
which means that we have entered the 3D-TF regime. If x ^ Ij we obtain 



1^^(3X)-. 



(12) 



This is the same result as is obtained from the stationary ID-GP equation using TF 
approximation. It means that we enter the ID mean-field regime. For the above local 
density approximation to be valid, Z ^ az must also be satisfied, which implies that 
X ^ 3^^^"^- This gives the boundary between the TF axial density profile and the 
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Figure 1. Phase diagram of the elongated condensate in the plane NX vs a±/as. The 
dashed lines indicate the crossover region between the ID mean-field regime and the 
ideal gas regime for A = 9.8 x 10^^ and A = 9.7 x 10^^, respectively 



Gaussian density profile (ideal ID gas). If the condition A^A <^ (os/o±)^ [11] is satisfied, 
the Tonks-Girardeau (TG) regime [23] can be reached. In this paper we will never enter 
this regime. 

Fig.l shows a plot of A^A vs a±/as along with several lines that schematically 
illustrate the phase that the gas is in. While the transition from one phase to another is 
gradual, these lines are a guide as to where the crossover is. As x decreases the gas will 
move from the 3D-TF regime to the ID mean-field regime. The red solid line on which 
A^A = a±/as, i.e. x = l? is marked to indicate the crossover region. As x decreases 
further the gas will enter the TG gas regime or the ideal gas regime. The ability to 
reach the ideal gas regime depends on the value of A and a±/as- The blue solid line, 
A^A = (cts/aj.)^, shows the crossover between the TG gas regime, and either the ID 
mean-field regime or the ideal gas regime. The ideal gas regime and the ID mean-field 
regime are divided by the line A^A = ^X^^'^a±/as. In Fig.l the location of this line is 
shown for two different values of A by the dashed lines. We can see that the 3D-TF 
and the TG gas regions move further far apart as a±/as becomes larger, leaving a larger 
parameter range for the ID mean-field (the ideal gas) description, and for certain values 
of a±/as and A, the ideal gas regime becomes inaccessible. The system will not enter 
the TG regime for the parameters looked at in this work {a±/as = 105.4, or 209.6 when 
A^A is greater than 10~^ as shown by the dotted lines in Fig.l). 

We consider an model, of which ®^Rb atoms are confined in a highly anisotropic 
harmonic trap. The trap frequency in the radial and axial direction is u}± = 27r x 91Hz 
and Uz = 2tt X 8.9Hz, respectively. The 3D s-wave scattering length is ag = 5.4nm. 
The system then undergoes a ID ideal gas, ID mean-field gas and 3D cigar-shaped gas 
transition with respect to an increasing number of atoms. 

The ground state of the condensate is obtained by the imaginary-time-evolution 
grid method for the time dependent effective-GP equation (5) and the 3D-GP equation 
(1). Theoretically the initial conditions we chose for the numerical simulation will only 
affect the efficiency of the calculation but not the final results. We found, however, 
that for the effective-GP equation (5), the choose of the initial chemical potential /i and 
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Figure 2. Relative differences of width (a) and peak density (b) of tlie axial condensate 
vs logioX calculated by the effective ID- and 3D-GP equation. 



the wave function during the transition between ID condensate and 3D cigar-shaped 
condensate will greatly affect the accuracy of the ground state properties and then 
the excitation frequencies. We compare our simulation results obtained by using the 
effective ID-GP equation (5) and the 3D-GP equation (1). As shown in Fig.2, the 
effective ID-GP equation gives the peak density and the width of the condensate with 
accuracy better than 2% in all regimes up to some critical value of x- Especially in the 
two limiting cases the results are nearly the same as those of the 3D-GP equation. The 
most interesting finding is that the variation of these two quantities is not smooth from 
the ID mean-field condensate to the 3D cigar-shaped condensate. In the crossover region 
there is a oscillation. We also find that there is a critical value of x where beyond this 
value the deviation of the results increases very quickly, which means that the effective 
ID-GP equation is not reliable anymore. This value is relevant to the aspect ratio of 
the trap. The larger the value of the aspect ratio the smaller the critical value. 



3. Bogoliubov equations 

In the frame of Bogoliubov theory, the gas is initially assumed to be in thermal 
equilibrium and the Bose field operator (j){z, t) can be expanded linearly as 

4>iz,t) = [M^) + 6^iz,t)]e~"^'/^ 

= 0000 + ^[uj{z)aj{t) - Vj{zya]{t)], 



(13) 



a'j, aj 



with 00 being the time independent groundstate amplitude of the condensate 
being the quasiparticle creation and annihilation operators, and where (Sq^o) = A^ is 
the number of atoms in the condensate. Uj and Vj are the qusiparticle amplitudes. 
Substituting this expression into the effective ID-GP equation, we can get the 
quasiparticle dynamical equation 



i?idt5(j){z, t) 



2m ' 



z, t) + Veff{z)6(j){z, t) + geffN<j)l6(j)^{z) (14) 
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where 

V{z)-^l + 2g,fJN\Mz)? ID 

Veff{z) = { (15) 

and 



V{z) - /i + Tiuj^Jl + 4a,iV|0oP + QeffNcpl effective-ID 



giD ID 

9eff = { (16) 

-^=^m effective-ID. 

To get these resuhs we have ignored all the higher order terms of 0{5^'^) and assumed 
that 5(j) + 5(f)^ is a very small quantity. The nonlinear term can be expanded as 

Jl + 4a,iV|0P^ Jl + 4a,iV|0oP+ , ^"-^'^^ =(^0 + ^0t). (17) 

/l + 4a,Ar|0o|2 



It is obvious that the ID result is nothing but the limiting case of the effective-ID result. 
It is then easy to show that Uj{z) and Vj{z) fulfill the Bogoliubov equations, 

/ A{z) -B{z)\(u,{z)\_ ( u,{z)\ 

[b*{z) -Aiz) )[v,{z) ) ""'[v.iz))' ^^'^ 

where the characteristic matrix elements are A(z) = —^d"^ + K//(-2^) and B{z) = 
geffN(j)l{z). By imposing the Rose commutation rules to the operators Oj and a^, we 
find that the quasiparticle amplitudes u and v must obey the normalization conditions 

dz[u*Uk - VjVk] = 6jk, / dz[u*Vk - ulvj] = 0. (19) 

In general, we can expand the mode functions u{z), v{z) and the condensate wave 
function (f){z) on a basis set ipn{z), 

00 00 00 

^i(^) = Y. (^n^V^n, Vj{z) = Y, fe^fVn, (^j{z) = ^ C^fVn, (20) 

n=0 n=0 n=0 

Practically, we can only involve a limited number of basis functions in our simulation. 
For an effective-ID BEC described by Eq.5, we denote the nonlinear interacting 
term as 

An,n' = (Jeff / dzj^ 4^^ 4^^^ 'Pn{z)ipi{z)(pi>{z)ipn'{z), (21) 

•' 1,1' 

and then Eq.(18) can be rewritten as 

where geff = geffKo-zh^z), Ej = Ej/hu^ X^y = {^n - A) <^n,n' + 2A„,„/, Yn,n' = ^n,n', 

and Ej is the eigenvalue of Bogoliubov equations. The value of En depends on the specific 
form of the basis function ipn- If we choose the simple-harmonic-oscillator (SHO) basis, 
which are the eigenstates of i^o = ~^|/2 + V^(5), then we get Ho'fn = {'n'+l/2)ipn = ^n'fn 
with n = 0, 1, . . . being the corresponding quantum number. Another choice is the 
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plane-wave (PW) basis, which is computationally the simplest (by virtue of the Fast 
Fourier transform method) in which En = ^^/2 + V^"^ with kn being the dimensionless 
momentum and V^" the Fourier part of \^ = 5^/2. 

From Eq.(22), we can get the condition for a nonzero solution is 

^^"^ (23) 



y-'Z \r-Z ^ 77i2 r 

A 



with / being a unit matrix. 

The quantum number n is defined as the number of nodes found in the axial wave 
function (pn{z). If there exists a solution of the normal-mode equations for positive Ej, 
then a negative solution also exists because if we replace Ej by —Ej and Uj by Wj, the 
Bogoliubov equations (18) remain unchanged. 

A formal solution of (18) at zero energy is -Eq = with uq{z) and vq{z) being 
proportional to the ground-state wave function of the condensate, but the norm of this 
solution is zero which is not normalizable in this way (19). The existence of a zero mode 
is consistent with Goldstone's theorem [24] since the nonzero average of the wavefunction 
breaks the U(l) global gauge symmetry of the Hamiltonian. 

To solve Eqs.(18) or (22) we need to get the ground state wave function of the 
condensate by solving the effective ID-GPE numerically. Then the coefficients a„, 6„ 
and the excitation energy, Ej, are calculated by matrix diagonalization providing a given 
basis set. 

The accuracy of the SHO method depends strongly on the number of basis states 
retained. The higher the energy modes we want to get, the more basis states we need to 
use. However, the largest number of basis states we can use are limited by the ability of 
solving the Hermit polynomials. So the number of excitation modes can be calculated 
is limited. Another limitation of this method is that to get accurate values for A, the 
integrand in Eq.(21) must be expressed as a polynomial. The basis functions V9„,(2;), for 
a harmonic trapping potential V{z), can be given by the dimensionless SHO functions 

1/2 



^n{z) 



-"'Hr.{~z). (24) 



By introducing another variable C, = z/\/2 and redefine the Hermit polynomial, Hn{z) 
the effective ID interacting term A„ „/ turns out to be a integral 



^(i) *(i) 



Kn' = T.^^^ J d(j9effe-^'H^iC)Hi{C)HiiC)H^iC). (25) 



w w' / J/- -■- 

1,1' ^2 

which can be evaluated by using a gaussian quadrature integration technique [25]. For 
the ID-GP equation, ^g// is just a constant which can be taken out of the integral, while 
for the effective ID-GP equation it is a function of the ground state wave function. Hence 
this method is not very accurate for the effective ID-GP equation (5). 

The PW expansion, related to the coordinate space, requires many more modes 
than the SHO basis. The size of the system and the mesh spacing are set when we 
calculate the ground-state wave function. The number of excitation modes we get is 
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the same as the number of mesh points, but only the first 20% are stable (accurate) 
when we change the size or number of mesh points. We will use the numerical results 
regarding to PW basis in our following discussions. 

4. Theoretical Analysis of collective excitations 

In a uniform 3D gas, the amplitudes u and v are plane waves and the resulting dispersion 
law for the elementary excitations takes the famous Bogoliubov form [26] 

1/2 



E{p) = hu(p) 



3Pp2 + (i: 

m \ 2m 



(26) 



where p is the momentum of the excitation and p is the density of the gas. For small 



momenta the spectrum takes the phonon-like form E{p) = cp, with c = Jgp/m being 
the sound velocity. At large momenta the quasiparticle behaves like free particle with 
energy E{p) ^ p^ /2m + gp. 

4.1. Variational analysis of the low-lying excitations in the low density domain 

For a ID-BEC the spectrum of the low energy excitations can be studied by following 
a time- dependent variational method introduced in Ref.[8] for the 3D case. 
Let's start with a dimensionless time- dependent GP equation 

-l^| + V{~z) + ~g^nN\4>\''^ 



idi(t){z,t) 



0(5, t). (27) 



2 

The problem of solving Eq. (27) can be restated as a variational problem corresponding 
to the minimization of the action related to the Lagrangian density 

^ = ^(M-0* - ~rdf^) + \d-A*di^ + y(5)|0|2 + ^A|0Y- (28) 

It is natural to choose a gaussian ansatz wave function 

-'i^+i2/3(i)+ii2^(t) 



0(5, t) = A{t)e~^^^^'''''"^'' ''" (29) 

in the very weak interaction regime, and the normalization is given by / \(f)\'^dz = 1. 
Here a,/?, 7 and a are all real parameters. The imaginary terms of the exponential 
are introduced to make the results reliable [8]. We insert Eq.(29) into Eq.(28) and 
calculate an effective Lagrangian, L, by integrating the Lagrangian density over the 
whole coordinate space L = /^ J^dz. Through a long but straightforward calculation, 
we get the equations of motion 

d?a + a = (30) 

d?a + a = — + ^. (31) 

The only one stable equilibrium point of Eqs.(30)-(31) is 

1 G 

tto = 0, o-Q 3 2=0 (32) 

'^0 "^0 
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E' 



2 



hu,, (34) 



where G = ^idA^/a/Stt. Expanding Eqs.(30) and (31) around the above equihbrium 
point (see Appendix ??), we can obtain the following frequencies for low-energy 
excitation modes: 

El = hu,, (33) 

, 3 G 1^/2 
+ 2^2^ + ^) 
where we have restored all the units. 

El corresponds to the dipole oscillation characterizing the motion of the centre of 
mass unaffected by the interatomic force. The oscillation frequency is just the frequency 
of the harmonic trap in the axial direction due to the harmonic confinement, which 
indicates that the motion of the centre of mass can be exactly decoupled from the 
internal degrees of freedom of the system. 

E2, the axial breathing mode, is the excitation energy corresponding to the width 
and the interaction strength of the condensate. The width of the condensate increases 
with the strength of the interaction between the particles. We note that Eq.(34) is only 
applicable in a restricted region where an ideal gas condensate with a Gaussian density 
profile is realized. In this region, E2 decreases with increasing G because the condensate 
expands rapidly in the axial direction. 

4-2. Analytical solutions in the hydrodynamic limit 

The analytic solutions of the linearized GP equation in the hydrodynamic limit for the 
low-energy axial modes were obtained in Refs.[27, 28]. The excitation frequency of the 
jth mode is given by 

Wj-iD = \/j(j + l)/2 ujz, (35) 

u}j-3DTF = yj(7+3)/4 u^, (36) 

for the ID mean- field and 3D-TF regime respectively, where j is a positive integer 
indicating the different excitation modes. These two equations are only valid if Uj <^ uj±. 
For noninteracting particles in a harmonic potential the ID form of the frequencies of 
the excitations are given by a; = juz- The dispersion relation of the normal modes 
of the condensate has changed significantly from the noninteracting behaviour, as a 
consequence of two-body interactions. Notice that in the case j = 1, the excitation 
frequency coincides with the oscillator frequency in the two limiting cases. According 
to the general considerations discussed in the previous section, this also holds for the 
crossover regime. The results of these low-energy modes, which can be directly excited 
by suitable modulation of the harmonic trap, have been confirmed by experiments with 
high precision[29, 30]. The fact that the frequency of the first excitation is not affected 
by the interatomic interactions, demonstrating that the lowest excitation in the system is 
the dipolar oscillation, offers a direct test on the numerical accuracy of our calculations. 
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Figure 3. The frequency of the three lowest excitation mode with respect to logwx- 
(curves from bottom: j = 1, 2, 3.) The arrows indicating the transition between the 
low density domain and the high density domain (see details in text). 



4-3. Collective excitations in the high density domain 

Following a sum-rule approach [4, 5, 31, 7], the collective frequency in the high density 
domain (from the ID mean-field regime to the 3D-TF regime) can be derived by using 
the formula [11] 

where {z^) = N~^ J_^ dzz'^lcpQ^z)]'^ is the expectation value of the square radius. This 
equation gives the same bounds to the frequency of the axial breathing mode as Eqs.(35) 
and (36). With Eqs.(7) and (10), we obtain 

3(7X1)3 -6x[(v^^)2 + 5] ' 
which means that the excitation frequency can be solely related to the axial half-length 
of the elongated condensate in the equilibrium configuration because x is also a function 
of Z (see Eq.(lO)). This provides a straightforward method to measure the excitation 
frequency experimentally. 



00 



5. Numerical results and discussion 

In Fig. 3, we show the Bogoliubov spectra of the three lowest excitations with respect 
to the parameter log^o X- The first excitation frequency is always Uz as predicted by 
the theory, while the frequencies of the second and the third excitation modes decrease 
gradually with x- The upper and lower bound of the axial breathing mode are ySuz and 
\/2.5ujz respectively, as calculated by Eqs.(35) and (36) when j = 2, while the frequency 
for the ideal gas limiting case is 2uz- The relative error of the frequency of the dipole 
mode is less than 0.15%, which means that our simulation results are with very high 
accuracy. Our simulation results also show a very good performance of recovering the 
same bounds of ID and 3D cigar-shaped condensate limiting cases for the second and the 
third excitation modes. However, the bounds of ID mean-field condensate, U2 = \/2^Uz 
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Figure 4. (a) The energy spectrum for the ID mean-field regime. The solid curve is 
the theoretical prediction from Eq.(35) and open symbols are numerical results (open 
squares, triangles and circles correspond to x = 0.47, 0.51, 0.93, respectively), (b) The 
energy spectrum for the 3D-TF regime. The solid curve is the theoretical prediction 
from Eq.(36) and the open symbols are the numerical results, (open squares, triangles 
and circles correspond to: x = 4.43 x 10^, 1.31 x 10^, 2.33 x 10^) 



and cua = -v/Ga;^ ..., are not completely recovered. Instead, there is a kink between the 
two domains (see arrows in Fig. 3). We shall explain this later. 

As shown in Fig.4, for the ID mean-field and the 3D-TF limiting cases, the 
frequencies of the very low-lying excitations with different interaction strengths are 
so close that they can be well approximated by a function of the excitation modes 
j (Eqs.(35) and (36)), as if they do not depend on the interatomic interaction. This 
differs from the uniform case where the dispersion relation, in the corresponding phonon 
regime, depends explicitly on the interaction through the velocity of sound. For a fixed 
value of X the accuracy of predictions (35) and (36) decreases as j increases, and in the 
ID mean-field regime, the frequency of the low-lying excitations is more sensitive to the 
variation of x than in the 3D-TF regime. 

In Fig. 5 the excitation spectra for different values of x ci-re plotted for one hundred 
excitation modes. The frequency of the excitations increases monotonically with j for 
a given x, while for a given j the frequency increases as x is reduced. The different 
spectra diverge as j increases. In the ID mean-field regime, this divergence decelerates 
with increasing j, whereas for 3D-TF regime it accelerates. 

In the low density domain (from the ideal gas regime to the ID mean-field regime) 
the theoretical prediction of the frequency of the axial breathing mode has been obtained 
as Eq.(34) by the variational analysis. In Fig. 6 we use this equation and the width of 
the condensate in the axial direction calculated by the effective ID-GP equation, ID-GP 
equation and 3D-TF approximation to compare the results with the results given by the 
Bogoliubov calculation. Both the results from the effective ID-GP equation and the 
ID-GP equation work well in this region. 

We use this equation and the axial half length of the condensate calculated through 
the effective ID-GP equation to obtain the frequency of the excitations in this region. 
Fig. 6 shows that our simulation is in good agreement with this prediction. 

As show in Fig.3, the Bogoliubov spectra recover the same bounds for the ideal 
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Figure 5. The energy spectrum vs the excitation modes, u^ — 2-k y. 8.9Hz, 
cji = 27r X 91Hz. (curves from top N=10, 100, 10^, lO", 10^ 10^.) Inset is a zoom 
plot for the very low-lying excitations. 



gas regime and the 3D-TF regime. However, there is a kink between the two domains 
(see arrows in Fig.3). The zoom plot of this region is in Fig. 6(a). As noted in Refs. 
[11, 12] there is a junction between the two domains originating from the finite value 
for a±/as- Keeping in mind that A <^ 1, then in order to have a perfect junction at 
{u2/ujz)'^ = 3, we should take the limit a±/as — )■ oo. This means that a well defined 
plateau corresponding to the ID mean-field regime can be fully reached from both sides. 
For aj_/as ~ 1 one would have a smooth transition between these two domains, without 
any plateau. Eq.(38) is valid in the high density domain, while Eq.(34) is valid only 
in the low density domain. We compare the numerical results with the theoretical 
predictions from Eqs.(38) and (34) for the two domains in Fig.6 and find they coincide 
very well, except when the system approaches the ID mean- field regime. The validity 
of the ID-GP equation is strongly related to the value of x- In Fig-6(a), we can see that 
the frequency of the excitation decreases dramatically after log^g X > 3, meaning that 
this method is not applicable any more in this region. 

6. conclusions 



In this paper we have calculated the Bogoliubov excitation spectra of a Bose gas in a 
very elongated trap. The system undergoes a ID ideal gas, ID mean-field gas and 3D 
cigar-shaped gas transition with respect to an increasing number of atoms or increasing 
the aspect ratio of the trap. In order to get the excitation spectrum for all these regimes, 
we take the effective ID-GP equation developed by Mateo et al. [2, 1] as a starting point. 
The Bogoliubov equations are solved by using the matrix diagonalization method in the 
plane wave basis. The results of our simulations are compared with those from the 
variational analysis in the low density domain, and the sum rule approach in the high 
density domain. We find that the Bogoliubov method fails to give the right spectrum 
in the transition region of these two domains. The plateau area is replaced by a kink 
with an irregular trajectory. We also find that there is a critical value of x where the 
effective ID-GP equation is not applicable any more. 
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iog,„); 



Figure 6. The frequency of the axial breathing mode, (a) The solid line is the 
theoretical prediction from Eq.(38). Open circles are numerical results, (b) The dashed 
line, dash dotted line and dotted line are the theoretical predictions from Eq.(34) 
by using the width of the condensate which is obtained from the ID-GP equation, 
effective ID-GP equation and 3D-TF approximation, respectively, a;^ = 27r x 8.9Hz, 
Wi = 27r X 91Hz. 
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